Cluster decay in very heavy nuclei in Relativistic Mean Field 
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Exotic cluster decay of very heavy nuclei has been studied in the microscopic Super- Asymmetric 
Fission Model. Relativistic Mean Field model with the force FSU Gold has been employed to obtain 
the densities of the cluster and the daughter nuclei. The microscopic nuclear interaction DDM3Y1, 
which has an exponential density dependence, and the Coulomb interaction have been used in the 
double folding model to obtain the potential between the cluster and the daughter. Half life values 
have been calculated in the WKB approximation and the spectroscopic factors have been extracted. 
The latter values are seen to have a simple dependence of the mass of the cluster as has been 
observed earlier. Predictions have been made for some possible decays. 
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• Decay of cluster of nucleons from very heavy nuclei was suggested by Sandulescu et ai. [1( in 1980. 
[T i | Subsequently, emissions of various clusters from very heavy nuclei leading to daughters around the 
magic nucleus 208 Pb have been observed^ experimentally. This new type of radioactivity is 
analogous to alpha decay where the decaying particle tunnels through the potential barrier. A lot 
of work has been done by various authors to calculate the half life for these exotic decay processes. 
Super-Asymmetric Fission Model (SAFM) has been used in [H and Preformed Cluster Model in 
[J] to describe the various observed decays. 

In the present work, we have used the microscopic SAFM. In most of the previous works in this 
model, the potentials that have been used were phenomenological in nature. Even in cases where 
the potential has been constructed from nuclear densities, the densities themselves were obtained 
using phenomenological prescriptions. However, microscopic densities obtained from mean field 
approaches may be expected to provide a better description of the densities and hence that of the 
' process of decay. With this in mind, we have chosen Relativistic Mean Field (RMF) to calculate 
the densities and study cluster decay in the present brief report. 

There has been a number of calculations for half life of cluster decay under the SAFM and we cite 
only a few recent ones. Basu Q has studied the nuclear cluster radioactivity also in the framework 
of a SAFM using a phenomenological density and the realistic M3Y interaction. Bhagwat and 
Gambhir Q have used densities from RMF calculations using the force NL30. In the present 
work, we have used another form of the Lagrangian density for RMF to study the structure of the 
ground state of the nuclei. We also have selected a form of the interaction with an exponential 
density dependence, DDM3Y1, different from that used in Ref. [f|. This particular form of the 
interaction is consistent in the sense that it can also explain the low energy scattering in the double 
folding approach. 

The interaction potential between the decaying cluster and the daughter nucleus has been ob- 
tained in the present work applying the double folding model by folding the densities of the cluster 
and the daughter nucleus using some suitable nuclear interaction along with the Coulomb interac- 
tion. As already mentioned, usually the densities are obtained from phenomenological description 
while in the present work, we have employed RMF to obtain the densities of the cluster and the 
daughter nucleus. 

RMF is now a standard tool in low energy nuclear structure, successfully reproducing various 
features such as ground state binding energy, deformation, radius, excited states, spin-orbit split- 
ting, neutron halo, etc. It is well known that in nuclei far away from the stability valley, the single 
particle level structure undergoes certain changes in which the spin-orbit splitting plays an impor- 
tant role. Being based on the Dirac Lagrangian density, RMF is particularly suited to investigate 
these nuclei because it naturally incorporates the spin degrees of freedom. In the case of exotic 
decays, certain neutron rich clusters such as 26 Ne, 30 Mg and 34 Si are emitted. We expect RMF 
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to successfully describe the nucleon densities in these clusters. There exist different variations of 
the Lagrangian density as well as a number of different parameterizations in RMF. Recently, a 
new Lagrangian density has been proposed [1] which involves self-coupling of the vector-isoscalar 
meson as well as coupling between the vector-isoscalar meson and the vector-isovector meson. The 
corresponding parameter set is called FSU Gold[8j. This Lagrangian density has earlier been em- 
ployed to obtain the proton nucleus interaction to successfully calculate the half life for proton 
radioactivity In this work also, we have employed this force. 

In the conventional RMF+BCS approach for even-even nuclei, the Euler-Lagrange equations are 
solved under the assumptions of classical meson fields, time reversal symmetry, no-sea contribution, 
etc. Pairing is introduced under the BCS approximation. Since the nuclear density as a function of 
distance is very important in our calculation, we have solved the equations in co-ordinate space. The 
strength of the zero range pairing force is taken as 300 MeV-fm for both protons and neutrons. 

The microscopic density dependent M3Y interaction (DDM3Y) may be obtained from a finite 
range nucleon nucleon interaction by introducing a density dependent factor. This class of in- 
teractions has been employed widely in the study of nucleon-nucleus as well as nucleus-nucleus 
scattering, calculation of proton radioactivity, etc. In this work, we have employed the interaction 
DDM3Y1 which has an exponential density dependence 

v(r, pi,p 2 , E) = G(l + aexp {-(3{ Pl + p 2 ))(l - 0.002E)u M3Y (r) (1) 

used in Ref. [l(| to study alpha- nucleus scattering. It uses the direct M3Y potential u M3Y (r) based 
on the G-matrix elements of the Reid[ll| NN potential and reproduces the saturation properties of 
cold nuclear matter. The weak energy dependence was introduced fl2j to reproduce the empirical 
energy dependence of the optical potential. The parameters used have the standard values viz. 
C = 0.2845, a = 3.6391 and (3 = 2.9605fm 2 . Here p\ and pi are the densities of the a-particle and 
the daughter nucleus, and E is the energy of the a-particle per nucleon in MeV. This interaction 
has been folded with the theoretical densities of cluster and the daughter nucleus in their ground 



states using the code DFPOT 13]. The assault frequency has been calculated from the decay 
energy following Gambhir et al 14J . The centrifugal barrier has been included in the cases where 
the I value for the decay is known to be non-zero. 

The decay probabilities have been calculated in the WKB approximation for penetration of the 
barrier by the cluster. No theoretical calculation can reproduce the Q-values for the decay very 
accurately. As the decay probability changes very rapidly with Q- value, we have taken the Q- values 
(and the decay energies) from experimental measurements. 

The spectroscopic factors have been calculated as the ratios of the calculated half lives to the 
experimentally observed half lives. This represents the preformation factor and may be considered 
as the overlap of the actual ground state configuration and the configuration representing the 
cluster coupled to the ground state of the daughter. Obviously it is expected to be much less than 
unity. The results for the spectroscopic factors are presented in TABLE [IJ Theoretical calculations 
have been performed for spectroscopic factors for some of these nuclei. Our results are comparable 
to those values. For example, our calculated value for S of 212 Po is 1.88 x 10~ 2 as compared to 
theoretical value 2.5 x 10 -2 deduced in (l5| . A value of 3.1 x 10~ 2 was obtained by Mohr 16] in a 
double folding model calculation using density from experimentally known charge distribution. 

It has been suggested[T3] that in the case of decay of heavy clusters, the spectroscopic factor 
may scale as 



S=(S a ) iA - 1 ^ 3 (2) 

where A is the mass of the heavy cluster and S a is the spectroscopic factor for the a-decay. Thus 
a plot of log 10 S against A should be a straight line. In Fig. [H we have plotted the negative of 
log 10 S for the decays where both the parent and the daughter are even-even nuclei against the 
mass number of the cluster and plotted a best fit line. One can see that the points fall nearly on a 
straight line with the S a value given by 1.93 x 10~ 2 . This is comparable to the value 1.61 x 10 -2 
obtained by Poenaru et aZfl8l]. 

We have also extended our study to decays where both the parent and the daughter nuclei have 
odd mass though the number of observed decays is rather small. The corresponding curve is shown 
in Fig. Here the S a value given by 1.35 x 10~ 2 . 

With such a good linear fit of the logarithm of spectroscopic factors with mass numbers, we have 
extended our scheme to calculate the half life of some other possible decays where unambiguous 
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lifetime measurements are not yet available possible or where there are possibilities of some decay 
taking place. Our results obtained with the fitted values of S from eqn. (1) are tabulated in 
TABLE |n] along with whatever experimental data are available. Except in the case of the decay 
of 233 U, the results are consistent with experimental observations. Even in the case of 233 U, the 
error is small. 

To summarize, we have studied exotic cluster decay of very heavy nuclei in the microscopic 
SAFM. The densities of the cluster and the daughter nuclei have been obtained from RMF. The 
microscopic interaction DDM3Y1 and the Coulomb interaction have been folded with these densi- 
ties to obtain the potential between the cluster and the daughter. The lifetime have been calculated 
in the WKB approximation and the spectroscopic factors have been extracted. The spectroscopic 
factors are seen to follow a simple rule as have been observed earlier. Predictions have been made 
for some possible decays. 

This work was carried out with financial assistance of the Board of Research in Nuclear Sciences, 
Department of Atomic Energy (Sanction No. 2005/37/7/BRNS). 
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TABLE I: Spectroscopic factor (S) of cluster decay obtained in the present calculation. 



Parent Cluster Q (MeV) log 1u T e3: (s) S 



Po 


4tj 

He 


8.95 


-6.52 


l.88x 10 


-2 


213d 

Po 


4tt 

He 


7.833 


-5.44 


1.67x 10 


-2 


214-t) 

Po 


4tt 

He 


7.833 


-3.78 


3.45x 10 


-2 


215 a j. 

At 


4tj 

He 


8.178 


-4.00 


1.31x10" 


-2 


221 ri 


14 

c 


31.317 


14.52 


1.50x 10" 


-8 


221 d 

Ha 


c 


on on/ 1 

32.396 


13.39 


1.55x 10 


-8 


222d 

Ra 


14^i 

c 


33.05 


11.00 


1.64x 10" 


-7 


223 d 

Ra 


c 


31.829 


15.20 


of .,m - 

2.85x 10 


-9 


224 d 

Ra 


c 


30.54 


15.92 


1.04x 10 


-7 


225 a 

Ac 


c 


OA A I-7T 

30.477 


17.34 


8.14x 10 


-8 


226d 

Ra 


C 


28.20 


21.34 


3.97x 10 


-8 


228^1, 

in 


o 


44.72 


20.72 


ot, . 1 n- 

8.37x10 


■11 


230 tt 
U 


22 t\t„ 

INe 


61.40 


19.57 


6.72x10 


■12 


230 rpv 

in 


24 at 

Ne 


57.571 


24.64 


1.87x10" 


■13 


231d 

Pa 


24i\t 

Ne 


60.417 


23.38 


3.13x10" 


■15 


232 tt 

u 


24h.t 

Ne 


62.31 


20.40 


9.77x10 


■14 


233 tt 


Ne 


60.486 


24.82 


1.47x 10 


■15 


234 tj 


24 Ne 


uo.ozu 


9^ 9^ 




■13 


233 tj 


25 Ne 


60.776 


24.82 


4.02x10" 


■16 


234 tj 


26 Ne 


59.466 


25.07 


1.67x10" 


■14 


234 tj 


28 Mg 


74.11 


25.74 


6.30x10" 


■17 


236p u 


28 Mg 


79.67 


21.67 


2.83x10" 


■17 


238p u 


28 Mg 


75.912 


25.70 


1.21x10" 


■16 


238p u 


32 Si 


91.19 


25.28 


2.34x10" 


■18 


238p u 


30 Mg 


91.19 


25.28 


2.34x10" 


■18 


242 Cm 


34 Si 


96.509 


23.15 


1.10x10" 


19 



TABLE II: Half life values of cluster decay obtained in the present calculation. 
Parent Cluster Q (MeV) logi T ej; (s) logioT^" 



^ 4 Th 


14 C 


32.929 




13.68 


226 Th 


14 C 


30.596 




18.28 


224 Th 


16q 


46.481 




15.47 


226 Th 


18 Q 


45.727 


> 16.8 


18.23 


232 Th 


24 Ne 


54.497 




29.96 


236TJ 


24 Ne 


55.945 




30.16 


232 Th 


26 Ne 


55.964 




28.57 


233tj 


28 Mg 


74.226 


> 27.6 


26.56 


237 Np 


30 Mg 


74.817 


> 27.6 


27.92 


240p u 


34 Si 


91.029 


> 25.5 


26.48 


241 Am 


34 Si 


93.926 


> 25.3 


26.25 




FIG. 1: Negative of logarithm of spectroscopic factors (S) as a function of cluster mass number A for 
even-even parents and daughters. 




FIG. 2: Negative of logarithm of spectroscopic factors (S) as a function of A for odd mass parents and 
daughters. 



